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Josephson effect in junctions of quasi one-dimensional triangular lattice superconductors is 
discussed, where the theoretical model corresponds to organic superconductors (TMTSF^PFg. 
We assume the quarter-filling electron band and p, d and / wave like pairing symmetries in 
organic superconductors. To realize the electronic structures in organic superconductors, we 
introduce the asymmetric hopping integral, (i') among second nearest lattice sites. At t' = 0, 
the Josephson current in the d wave symmetry saturates in low temperatures, whereas those 
in the p and the / wave symmetries show the low-temperature anomaly due to the zero-energy 
state at the junction interfaces. The low-temperature anomaly appears even in the d wave 
symmetry in the presence of t' , whereas the anomaly is suppressed in the / wave symmetry. 
The shape of the Fermi surface is an important factor for the formation of the ZES in the 
quarter-filling electron systems. 
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1. Introduction 

Quasi one-dimensional (Q1D) superconductors 
(TMTSF) 2 X (X=PF 6 , C10 4 , etc.) 1,2) have recently 
attracted much attention as a possible spin-triplet su- 
perconductor. A spin-triplet pairing has been suggested 
from an observation of large H C 2 3 ^ and unchanged 
Knight shift across T c . 4 ' Unconventional pairing with 
nodes on the Fermi surface has been suggested from 
an NMR measurement, 5 ' while a thermal conductivity 
measurement has reported the absence of nodes on the 
Fermi surface. 6 ' In addition, a recent experiment showed 
the zero-bias conductance peak (ZBCP) in tunnel 
spectra. 7 ' Theoretically, a p wave paring in which nodes 
of the pair potential can be made to avoid intersecting 
the Q1D Fermi surface has been proposed in an early 
stage. 8-10 ' On the other hand, a spin-singlet d wave pair- 
ing mediated by spin fluctuations has been proposed by 
several authors. n ~ 13 ' This is because superconductivity 
lies right next to the Ikp spin density wave (SDW) phase 
in the pressure-temperature phase diagram. Moreover, 
one of the present authors has recently proposed that 14 ' 
a triplet / wave pairing may dominate over the d and 
the p wave in (TMTSF^PFg due to a combination of 
a Q1D Fermi surface, coexistence of the 2kp SDW and 
the 2kp charge density wave (CDW) suggested from 
a diffuse X-ray scattering, 15 ' 16 ' and an anisotropy in 
the spin fluctuations. Thus the situation is not settled 
either experimentally or theoretically. In recent papers, 
we proposed the tunneling conductance experiment to 
determine which one of pairing symmetries is realized 
in (TMTSF) 2 X 17 ' 18) because the tunnel spectra are 
sensitive to the internal degree of freedom of Cooper 
pairs. 

A great variety in the Andreev reflection 19 ' is one of 



interesting features in superconductors with unconven- 
tional pairing symmetries. The pair potential of a trans- 
mitted quasiparticle in the electron branch differs from 
that in the hole branch because of unconventional pair- 
ing symmetries. When the two pair potentials have op- 
posite signs to each other, constructive interference ef- 
fects of a quasiparticle near the surface of supercon- 
ductors lead to the zero-energy state (ZES). 20-25 ' The 
ZES is observed as the ZBCP in tunneling conductance 
spectra. Since the high-T c superconductors have the d 
wave pairing symmetry, 26 ~ 31 ' the ZBCP was found in a 
number of experiments. 23, 24 > 32 ~ 44 ) So far a considerable 
number of theoretical studies have been made on the 
ZES itself 20 ' 21 ' 25 ' and related phenomena of transport 
properties in both spin-singlet 45-61 ' and spin-triplet 62-68 ' 
unconventional superconductor junctions. Effects of fer- 
romagnets attaching to superconductors, 69-76 ' those of 
broken time-reversal symmetry states (BTRSS), 77-89 ' 
those of magnetic fields 90-93 ' and those of the random 
potentials 94 ~ 105 ' on the ZBCP are hot topics in recent 
studies. Since the ZES appears just on the Fermi energy, 
it drastically affects transport properties through the in- 
terface of unconventional superconductor junctions. The 
low-temperature anomaly of the Josephson current be- 
tween the two unconventional superconductors is ex- 
plained in terms of the resonant tunneling of Cooper 
pairs through the ZES. 66 ' 103 ' 106-111 ' The Josephson ef- 
fect in high-T c superconductors has been reported in a 
number of papers. 112-122 ' 

In theories, the Andreev reflection in unconventional 
superconductors has been studied based on the free elec- 
tron model in which the quadratic dispersion relation 
and the isotropic Fermi surface are assumed. Although 
the Fermi surface in real unconventional superconductors 
are not isotropic at al, the free electron model well ex- 
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plains characteristic behaviors of the tunnel conductance 
and the Josephson current. This may be because such 
transport properties are sensitive to the pairing symme- 
tries of superconductors. However recent studies show 
several exceptions. 18 ' 60 ' For instance in high-T c super- 
conductor junctions, the free electron model does not 
explain the Josephson current calculated on the two- 
dimensional tight-binding model for some cases. 60 '' The 
low-temperature anomaly of the Josephson current is ex- 
pected in the free electron model when the a axis of the 
high-T c superconductors oriented from the junction in- 
terface normal. The anomaly is washed out by the Friedel 
oscillations of the wave function in the lattice simulation 
for specific orientation angles. In addition, the density of 
states at the surface of the superconductors is sensitive to 
the shape of the Fermi surface in Q1D organic supercon- 
ductors (TMTSF)2X when triangular lattice structures 
are taken into account. In such situation, the Joseph- 
son current is also expected to be sensitive to the pair- 
ing symmetries and the shape of the Fermi surface. In 
this paper, we discuss the direct-current Josephson ef- 
fect motivated by the surface density of states in Q1D 
superconductors. 17, 18, 123 ' So far, a Josephson effect in 
organic superconductors has been reported in the theo- 
retical paper. 124 ) The effects of electronic structures in 
QID lattice, however, are not taken into account. 

This paper is organized as follows. In Sec. 2, we de- 
scribed the Josephson junctions of organic superconduc- 
tors by the Bogoliubov-de Gennes equation on the two- 
dimensional lattice. The Josephson current is discussed 
for p, d and / wave symmetries in Sec, 3. In Sec. 4, we 
discuss the calculated results. We summarize this paper 
in Sec. 5. 

2. Model and method 

Let us consider QID superconductor / insulator / su- 
perconductor (SIS) junctions as shown in Fig. 1, where 
r = jx + my labels a lattic site, where x and y are unit 
vectors in the x and the y directions, respectively. The 
two superconductors (i.e., — oo < j < and L + 1 < j < 
oo) are separated by the insulator (i.e., 1 < j < L). We 
assume the periodic boundary condition in the y direc- 
tion and the number of lattice sites in the y direction is 
M . The junctions are described by the mean-field Hamil- 
tonian 
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(c r (T ) is the creation (annihilation) operator 
of an electron at r with spin a = ( | or j), &q is the 
2x2 unit matrix representing the spin space, frj with 




A-type (b) Mirror B-type 

i ®!7Qri®fO-® :<I>tQd3b<>!<I>: 

/?, d wave /wave 

Fig. 1. A schematic figure of the SIS junction of the organic su- 
perconductor. 



j = 1 ~ 3 are the Pauli matrices and S in the summa- 
tion denotes the superconductors. The Fermi energy in 
the superconductor is [i r = fis for j < and j > L + 1, 
and that in the insulator is \i r = hn for 1 < j < L. In 
superconductors, we assume e r = and t a and tb are the 
hopping integral in the x and y directions, respectively. 
We also introduce asymmetric hopping integral (i') be- 
tween the second nearest neighbors as shown in Fig. 1 
to realize electronic structures in (TMTSF^PFg. In the 
insulator, e r — Vb for 1 < j < L denotes the barrier 
potential and the hopping integrals in the two directions 
are equal to t a - The pair potential in p, d and / wave 
symmetries are defined by, 
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where (pi = ifL or ipn is the macroscopic phase of super- 
conductors and e.j with j = 1, 2 and 3 are unit vectors 
in the spin space. A schematic picture of Cooper pairs 
are shown in Fig. 1. The Hamiltonian is diagonalized by 
the Bogoliubov transformation, 
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where 7^ (7a, o-) is creation (annihilation) operator of 
a Bogoliubov quasiparticlc. In Eq. (8), u\ and v\ are 
the wavefunction of a quasiparticle which satisfy the 
Bogoliubov-de Gennes (BdG) equation. 125 ' 
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The eigenvalue E\ is independent of spin channels be- 
cause we consider unitary states in superconductors. In 
what follows, we briefly discuss the method to calculate 
the Josephson current for the d wave symmetry. The ap- 
plication to the p and the / wave symmetries is straight- 
forward. In the case of the d wave symmetry, the BdG 
equation in Eq. (10) is decoupled to two equations, 
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where (wy), for example, represents an element of u in 
Eq. (8) and [i^ijfn] obeys essentially the same equa- 
tion. In the following, we suppress (1,1) of («n) and 
(2, 1) of («2i)- In this way, the wave function 

/ u x (jx + ly) \ 



*A(j) = 



(12) 



u\(jx + My) 
v\(jx + ly) 

V v x (jx + My) j 

satisfies the BdG equation. For j < —2, for instance, the 
BdG equation reads 



A 







2 V e-^ L l 
' t N (-) 



+ 



+ 







6 

6 

-Tjv(- 



-Exl + Es 

6 -Exi - E s 

f N (+) 6 
6 -f N (+) 



*a(.? + 2) 

*a(j) 

- 1) 



A 



Tn(+) = 







e l<PLl 
6 



*a(7 - 2) = 0, (13) 



/ -t a ~t' 
-t„ -if 



\ -t' 

/ -t a ••• 
-if -ta 



Tn(-) 



\ 




-ta -t' 
-t a ) 

-if \ 





(14) 



••• - 

\ ••• 

/ -Ms -h 

-tb -/is -*6 



if -ta 



(15) 



^s = 



-a -t a J 
-tb \ 





o 

V -t b 



-h -^s 



-tb 



(16) 



where 1 and are the M x M unit matrix and the zero 
matrix, respectively. To solve the BdG equation, we ap- 
ply the recursive Green function method 126-128 ^ and cal- 
culate the Matsubara Green function in a matrix form 

GWO',/) =^* A (j)[^ rt -^A]- 1 *i(j'), (17) 



where w n = (2n + 1)ttT is the Matsubara frequency 
and T is a temperature. Throughout this paper, we 
use the units of h = ks = 1, where fee is the Boltz- 
mann constant. The Josephson current in the insulator 
(1 < j < L), is given by 127 " 129 ) 
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We note that J(j) is independent of j when we consider 
the direct-current Josephson effect. 

3. Josephson current 

The low-temperature anomaly of the Josephson 
current is a typical phenomenon in the quantum 
transport between two unconventional superconduc- 
tors. 66, 103 ' 106_111 ) In this section, we discuss effects of 
the asymmetric second nearest neighbor hopping (t') on 
the ZES at the interface and on the Josephson current. 
For finite t' , it is possible to consider two types of SIS 
junctions as shown in Fig. 1 (a) and (b). The parallel 
junction consists of two A-type superconductors, whereas 
the mirror-type junction consists of a A-type and a B- 
type superconductors. Generally speaking, the Josephson 
effect in the two junctions are not identical to each other. 

In what follows, we choose parameters as tb = 0.1t a , 
fi s = -1.4099i o , /i n = -2.0t a , M = 20, L = 4, and 
Vb = 2.0t a . The electron density is fixed at the quarter- 
filling. The amplitude of the pair potential at the zero 
temperature is Ao = 0.1t a and the dependence of the 
pair potential on temperatures is described by the BCS 
theory. 

3.1 p wave symmetry 

In Fig. 2, we show the Josephson current in the par- 
allel junctions with the p wave symmetry. The vertical 
axis is normalized by 7rA /2ei?]v, where Rn is the nor- 
mal resistance of junctions. We note in this vertical scale 
that the Josephson current of the s wave symmetry with 
the isotropic Fermi surface is close to unity in the limit 
of the zero temperature. 130 - ) In real materials, the second 
nearest neighbor hopping should be —t'/t a = 0.08. We 
also show results for —t'/t a = and 0.16 for compari- 
son. In Fig. 2 (a), the maximum value of the Josephson 
current (J m ax) is plotted as a function of temperatures, 
where (J ma x) is estimated from the Josephson current 
as a function of ip = tpL — <Pr- The Josephson current 
increases with decreasing temperature irrespective of t' 
and does not saturate even in low temperatures. This 
behavior is called as the low-temperature anomaly of the 
Josephson current and is owing to the ZES forming at 
the junction interface. When the d vector has only one 
component in the free electron model, a condition for 
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Fig. 2. Joscphson current for the p wave symmetry is shown, 
where t b = 0.1t a , MS = -1.4099t o , M = 20, L = 4, and 
Vb = 2t a . In (a), the maximum amplitude of the Josephson 
current (Jmax) is plotted as a function of temperatures. In (b), 
Jmax at T = 0.017c are shown as a function of —t'/t a . The 
current-phase relation is calculated at T = 0.01T C in (c). 



appearance of the ZES is given by 

d(k a ,h) ■ d(-k a ,k b ) < 0, 



(19) 



where d(k ai k b ) is the Fourier component of d r - r >, and 
k a and k b are the wavenumber in the x and y directions, 
respectively. The two pair potentials in Eq. (19) corre- 
sponds to the pair potentials in the electron and the hole 
branch of a quasiparticle. Since the translational invari- 
ance in the y holds, kb is conserved in the transmission 
and the reflection of a quasiparticle at the junction in- 
terfaces. In the p wave symmetry, the pair potential in 
Eq. (5) results in 

d {p \k a ,k b ) = Asin(2fc a )e 3 . (20) 

Since Eq. (20) is an odd function of k a , it satisfies 
Eq. (19) for all the Fermi surface as shown in Fig. 3, 
where we draw the Fermi surface for t' = in (a), 
t' = — 0.08i a in the A-type superconductor in (b), and 
t' = — 0.08t a in the B-type superconductor in (c) with 
the solid line. The pair potential in the p wave symme- 
try has a node line at k a = in both the free electron 
model and the lattice model. In the lattice model, how- 
ever, Eq. (20) has additional node lines at k a = ±0.57r be- 
cause the pairing interaction works between two electrons 
on the second nearest neighbor sites at the quarter-filling 
as shown in Fig. 1. We define the additional node lines 
as the node lines appear in the pair potentials because 
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Fig. 3. The Fermi surface of TMTSF for t' = is shown in (a). 
Those for t' = — 0.08t a in the A-type and the B-typc supercon- 
ductors arc shown in (b) and (c), respectively. The pair potential 
in the p wave symmetry is given by A sin 2k a which changes the 
sign at k a = as indicated by + and - in the figures. 



of the quarter-filled electron band on the Q1D lattice. 

The asymmetric hopping (£'), modulates the shape of 
the Fermi surface as shown in Fig. 3 (b) and (c). When 
the Fermi surface looses a symmetry with respect to k a — 
as in (b) and (c), the condition for the ZES should be 
rewrote as 

d(k a ,k b )-d(-k' a ,k b )<0, (21) 

where k a and — k' a are the wave numbers of the Fermi 
surface for fixed k b . In (b), the two wavenumbers are in- 
dicated by arrows on the Fermi surface for k b = tt/2. In 
the following, we call Eq. (21) as the zero-energy condi- 
tion (ZEC). Even in the presence oft', the ZEC is always 
satisfied because the additional node lines are far from 
the Fermi surface. The amplitude of the Josephson cur- 
rent at T — 0.01T C decreases with increasing —t'/t a as 
show in Fig. 2 (b). This is mainly because Rn decreases 
with increasing —t'/t a . In our calculation, Rn in units 
of h/e 2 are 1.04, 0.83 and 0.63 for -t'/t a = 0, 0.08 and 
0.16, respectively. In Fig. 2 (c), the current-phase relation 
is calculated at T = 0.01T C . The current-phase relation 
in a low temperature deviates from the sinusoidal func- 
tion because the resonant tunneling through the ZES en- 
hances the multiple Andreev reflection between the two 
superconductors. We note that the results in the mirror- 
type junction are identical to those in Fig. 2 (a)-(c). In 
the p wave symmetry, the Josephson current are not qual- 
itatively changed by introducing t' because the ZES at 
the interface governs the characteristic behavior of the 
Josephson current. 

3.2 d wave symmetry 

In Fig. 4, we show the Josephson current in the parallel 
junctions with the d wave symmetry. 

In (a), Jmax is plotted as a function of temperatures. 
In the absence of i', the Josephson current saturates in 
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Fig. 4. Joscphson current in the parallel junction with the d wave 
symmetry is shown. In (a), the Josephson current is plotted as 
a function of temperatures. In (b), amplitudes of the Joscphson 
current at T = 0.01T C are shown as a function of —t'/ta- The 
current-phase relation is calculated at T = 0.01T C in (c). 
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Fig. 5. Joscphson current for the d wave symmetry in the mirror- 
type junction. 
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Fig. 6. The Fermi surface of TMTSF for t' = is shown in (a). 
Those for t' = — 0.08t a in the A-type and the B-type supercon- 
ductors are shown in (b) and (c), respectively. The pair potential 
in the d wave symmetry is given by A cos 2k a which changes its 
sign at k a = ±0.25-7r as indicated by + and - in the figures. The 
Fermi surface for t' = — 0.Wt a is shown for A-type and B-type 
superconductors in (d) and (e), respectively. The Fermi surface 
shown with the solid line satisfy the condition for the ZES in 
Eq. (22). 



low temperatures as that in the s wave junctions. On 
the other hand, the results for —t'/t a = 0.08 and 0.16 
do not saturate even in low temperatures. In particular, 
the Josephson current for — t'/t a =0.16 rapidly increases 
with decreasing temperatures. In the case of spin-singlet 
superconductors, a condition for appearance of the ZES 
is given by 

d(k a ,k b )d(-k' a ,k b ) <0, (22) 

where d(k a , kb) is the Fourier component of d r - r > . In the 
d wave symmetry, the pair potential in Eq. (6) results in 

dW(fc ,fc 6 ) = Acos(2fc ), (23) 

and is the even function of k a . The additional node lines 
in the lattice model are k a = ±0.257r and ±0.757r. In 
the absence of t', Eq. (23) does not satisfy the ZEC in 
Eq. (22) for all the Fermi surface as shown in Fig. 6, 
where we show the Fermi surface for t' = in (a). The 
Fermi surface indicated by the broken line does not sat- 
isfy the ZEC. In the presence of t', however, the ZEC can 
be satisfied for some wave numbers on the Fermi surface 
as plotted with the solid line in (b) for t' = — 0.08£ a . The 
sign of the pair potential changes frequently on the Fermi 
surface because t' modifies the shape of the Fermi sur- 
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face lies along the additional node lines at k a = ±0.257r. 
As a result, the ZEC in Eq. (22) is satisfied for some 
wavenumbers, which leads to the formation of the ZES. 
This implies an importance of lattice structures on the 
quantum transport. The current-phase relation deviates 
from the sinusoidal relation as shown in the results for 
—f/t a = 0.16 in Fig. 4 (c) because most of the Fermi 
surface satisfy the ZEC as shown in Fig. 6 (d). We note 
that the Fermi surface around kb — and ±ir are still 
out of the ZEC for -t' /t a = 0.16. The deviation for 
—t'/t a = 0.08 from the sinusoidal relation is smaller than 
that for —t'/ta = 0.16. A quasiparticle incident perpen- 
dicular to the junction interface mainly contributes to the 
Josephson current. The Fermi surface at kb = 0, however, 
does not satisfy the ZEC as shown in Fig. 6 (b). Thus 
the contribution of the resonant tunneling via the ZES 
is small for —t'/t a = 0.08 even though the temperature 
dependence tends to be anomalous in low temperatures. 

The Josephson effect in the mirror-type junctions is 
different from that in the parallel junctions because of 
the lattice structures. In Fig. 5, we show the Joseph- 
son current for the d wave symmetry in the mirror- 
type junctions. In (a), J max for —t'/t a = 0.08 first in- 
creases with decreasing temperatures then decreases for 
T < 0.05T C . Such non monotonic temperature depen- 
dence has been also reported in the high-T c superconduc- 
tor Josephson junctions. For —t'/t a = 0.16, the Joseph- 
son current changes its sign and the amplitude increases 
rapidly with decreasing temperatures. The differences 
between the parallel and the mirror-type junctions can 
be understood in terms of the relative sign of the pair 
potentials in the two superconductors. In the presence 
of the time-reversal symmetry, the Josephson current is 
decomposed into a series of 

oo 

J =^J n sin(nip). (24) 

n=l 

It was shown that the J\ is roughly given by 111 ) 

Ji = ^2d R (k' a ,kb)d L (k a ,kb)F 1 (k b ), (25) 

k b 

where F\ is the positive function of kb, dn(k' a ,kb) and 
dL(k a ,kb) are the pair potential on the Fermi surface in 
the right and the left superconductors, respectively. Thus 
the product of du(k' a ,kb)dL(k aj kb) determines the sign 
of the Josephson current proportional to sin p. The par- 
allel junction consists of two A-type superconductors as 
shown in Fig. 1. Therefore d^(k a ,k b ) and di,(k a ,kb) are 
identical to each other, which results in the positive sign 
of J\. On the other hand in the mirror-type junctions, 
di(fco) in the A-type superconductor shown in Fig. 6 (b) 
has the opposite sign to dn(k' ) in the B-type supercon- 
ductor in (c). It is easily shown that dn(k' al kb)di,(ka, kb) 
is negative when the Fermi surface satisfies Eq. (22) in 
the mirror-type junctions. In high temperatures, a quasi- 
particle on the Fermi surface shown with the broken lines 
in Figs. 6 (b) and (c) dominates the Josephson current 
and Ji is positive. This is because a quasiparticle inci- 
dent perpendicular to the junction (kb — 0) mainly con- 
tributes to the Josephson current and the contribution of 
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the ZES is negligible. In low temperatures, the resonant 
transmission via the ZES also contributes to the Joseph- 
son current. As a consequence, Ji has a non monotonic 
temperature dependence as shown in Fig. 5 (a) because 
the sign of the Josephson current via the ZES is nega- 
tive. For —t'/ta — 0.16, the Josephson current is negative 
for low temperatures because the resonant transmission 
via the ZES dominates the Josephson current. It is also 
shown that the Josephson current proportional to sin 2ip 
is given by 

J 2 = -Y,iMka,kb)d L (k' a ,k b )} 2 F 2 (k b ), (26) 

where F 2 is the positive function of kb- The sign of J 2 
is always negative irrespective of the pairing symmetries 
of the two superconductors. Thus the current-phase re- 
lation for —t'/ta = 0.16 in Fig. 4 (c) takes its maximum 
at ip > 0.57T in the parallel junctions. While that in the 
mirror-type junction takes its minimum at ip < 0.57T in 
the mirror-j unction in Fig. 5 (c). 

3.3 f wave symmetry 

In Fig. 7, we show the Josephson current for the / 
wave symmetry in the parallel junctions. 

In (a), the maximum amplitude of the Josephson cur- 
rent is plotted as a function of temperatures for —t'/t a 
— 0, 0.08 and 0.16. In the absence of t' , the Joseph- 
son current show the low-temperature anomaly because 
the Fourier component of the pair potential is given by 
Asin4fc a which satisfies Eq.(21) for all the Fermi surface 
as shown in Fig. 9 (a). The additional node lines are at 
k a = ±0.257r, ±0.57r and ±0.757r in the lattice model. 
The anomalous behavior tends to disappear for finite 
t'. In Fig. 9 (b), the Fermi surface which do not satisfy 
Eq.(21) is indicated by the broken line for t' = — 0.08t a . 
For t' — — 0.16i a , most wave numbers of the Fermi sur- 
face do not satisfy Eq.(21) as shown in Fig. 9 (d). In 
contrast to the d wave junctions, the ZES is suppressed 
by introducing t' in the / wave junctions. The phase- 
current relation for —t'/t a = and 0.08 in Fig. 7 (c) 
apparently deviates from the sinusoidal function because 
of the multiple Andreev reflection via the ZES. The re- 
sults for —t'/ta = 0.16 slightly deviate from shop, but 
degree of deviation is smaller than those for —t'/t a = 
and 0.08. 

In Fig. 8, we show the Josephson current for the / 
wave symmetry in the mirror- type junctions. In (a), the 
Josephson current for —t'/t a = and 0.08 are essentially 
the same as those in the parallel junctions. In the current- 
phase relation for —t'/t a = 0.16 in (c), the Josephson 
current becomes zero around p = 0.757T. In this case, 
J\ is negative because dn(k' a , kb) ■ d^(k a , kb) is negative 
for almost all the Fermi surface. The amplitude of J 2 is 
not negligible because the resonant transmission through 
the ZES for a quasiparticle incident perpendicular to the 
interface (kb = 0) enhances J 2 . In a rough estimation, we 
find |J 2 | ~0.7|Ji|. 

4. Discussion 

The formation of the ZES is a universal phenomenon at 
the surface of unconventional superconductors. 25 ^ It may 
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Fig. 7. Joscphson current for / wave symmetry in the parallel 
junctions. In (a), the maximum amplitude of the Josephson cur- 
rent is plotted as a function of temperatures. In (b), amplitudes 
of the Josephson current at T = 0.01T C are shown as a function 
of —t'/ta- The current-phase relation is calculated at T = 0.01T C 
in (c). 
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Fig. 9. The Fermi surface of TMTSF for t' = is shown in (a). 
Those for t' = — 0.08t a in the A-type and the B-type junctions 
are shown in (b) and (c), respectively. The Fermi surface for 
t' = — 0.16i a is shown in (d) and (e). The pair potential in the / 
wave symmetry is given by A sin Ak a which changes the sign at 
k a = and ±0.25-7r as indicated by + and - in the figures. The 
Fermi surface shown with the solid line satisfies the condition 
for the ZES. The broken line indicates the wave numbers of the 
Fermi surface which do not satisfy the condition of the ZES. 



be possible to analyze the pairing symmetry of unconven- 
tional superconductors from the anisotropy in the tun- 
neling conductance and that in the Josephson current. 
Actually we have investigated these transport properties 
in both the free electron model and the lattice model. As 
far as we study, there are no drastic differences between 
the two theoretical models in the case of high-T c super- 
conductors. The formation of the ZES has been believed 
to be insensitive to the electronic structures such as the 
shape of the Fermi surface. The calculated results of this 
paper indicate that this statement is justified when the 
Fermi surface is far from the additional node lines pe- 
culiar to the lattice structures. We have shown that the 
shape of the Fermi surface strongly affects the formation 
of the ZES in Q1D organic superconductors. The Q1D 
quarter-filling band and the asymmetry of the Fermi sur- 
face with respect to k a are important factors for the ZES 
to be sensitive to the shape of the Fermi surface. In the 
quarter-filled band, the pairing interaction tends to work 
between two electrons at the second nearest neighbor 
sites in the x direction as shown in Fig. 1 because of 
the long range Coulomb repulsion. In particular, the / 
wave superconductivity requires the pairing correlation 
between the 4th nearest neighbor sites in the x direc- 
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tion. As a consequence, the pair potential changes sign 
at additional node lines in the first Brillouin zone such 
as k a — ±0.57r and ±0.257r as shown in Figs. 3, 6 and 
9. The asymmetric hopping (t') which characterizes the 
triangular lattice structure causes the asymmetry of the 
Fermi surface with respect to k a . The pair potential of- 
ten changes its sign along the asymmetric Fermi surface 
because the Fermi surface in Q1D quarter-filled band lies 
along the additional node lines at k a — ±0.257r. In the p 
wave symmetry characteristic behaviors of the Joseph- 
son current are insensitive to t' because the Fermi surface 
is far from the additional node line at k a — ±0.57r. The 
sign change of the pair potential around the additional 
node lines at k a = ±0.257r is very important for the for- 
mation of the ZES in the d and / wave symmetries. In 
the absence of t', it is clear that the ZES is formed in 
the / wave symmetry and that the ZES is not formed 
in the d wave symmetry. This is because the pair poten- 
tial in the / wave (d wave) symmetry is an odd (even) 
function of k a . Roughly speaking, t' assists the formation 
of the ZES in the d wave symmetry and suppresses the 
ZES in the / wave symmetry. As a results, the Joseph- 
son current exhibit the various temperature dependences 
and the phase current relations depending on the pair- 
ing symmetries, the degree of asymmetry and the types 
of junctions. 

5. Conclusion 

We have studied the Josephson current in quasi one- 
dimensional unconventional superconductors with trian- 
gular lattice structures. The theoretical model describes 
the organic superconductors such as (TMTSF) 2 X. In 
the calculation, we assume p, d, and / wave like pair- 
ing symmetries in superconductors. The pair potentials 
have the additional node lines because the pairing inter- 
action works between the second or 4th nearest neigh- 
bor sites in the current direction at the quarter-filled 
electron band. The triangular lattice structure is charac- 
terized by the asymmetric second nearest neighbor hop- 
ping (t') introduced on the square tight-binding lattice. 
The zero-energy states is sensitive to t' in the d and 
/ wave symmetries because the Fermi surface lies just 
along the additional node lines at k a = ±0.257r. It is 
possible to consider the two types of junctions, (i.e., par- 
allel and mirror), because of the triangular lattice struc- 
tures. The Josephson effect in parallel junctions is qual- 
itatively different that in the mirror-type junctions. We 
show that the Josephson current has various tempera- 
ture dependences and current-phase relations depending 
on the pairing symmetries of superconductors, the shape 
of the Fermi surface and the types of the junctions. 
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